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SUMMARY 


A method is presented for calculating the aerodynamic lifts and 
moments experienced by a cascade or two-dimensional approximation to a 
compressor or turbine blade row in subsonic flow under harmonic oscil- 
lation. Arbitrary stagger and interblade phase-lag angles are permitted. 
The most significant features of the method stem from the utilization of 
Fourier transforms of blade pressure- jump functions. This permits expres 
sion in closed form of the kernel function appearing in the integral equa 
tion relating upwash to the pressure-jump transform. So-called resonance 
phenomena, first discovered in connection with the subsonic wall inter- 
ference problem, are shown to occur in the case of cascade oscillation 
even in the presence of stagger and arbitrary interblade phase-lag angles 
The resonance points are shown to be related to the poles of the kernel 
function for the Fourier transform formulation of the problem. Numerical 
techniques are developed for the direct solution of the transform rather 
than of the pressure jump proper, and it is shown that lift and moment 
may be easily expressed in terms of the pressure-jump transform without 
requiring any inversion. The method is applied to computation of a zero- 
stagger cascade in antiphase motion (corresponding to the tunnel-wall 
interference problem) in order to permit comparison of results with those 
obtained by solution of the more conventional integral equation based 
upon the Hankel function series type of kernel. Excellent agreement is 
obtained in all cases. 


INTRODUCTION 


The aerodynamic forces and moments experienced by compressor blades 
undergoing oscillatory motion or in the presence of certain classes of 
oscillatory inflow are of interest to the aircraft -jet-engine industry 
for several reasons. First, while the occurrence of classical flutter 
in compressor blades has not been a matter of concern up to this time, 
it is conceivable that current design trends may lead to configurations 
and flow conditions in which the avoidance of classical flutter must 
become a consideration. Second, the value of aerodynamic damping is of 
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interest quite apart from the flutter problem proper, and the existence 
in compressible subsonic flow of certain so-called resonance phenomena 
can have a major influence upon this source of damping. Third, the 
aerodynamic behavior of blade .rows under the influence of certain types 
of inflow disturbance or nonuniformity may be calculated by the same 
methods utilized in the computation of aerodynamic forces and moments 
experienced by the blades undergoing oscillations. Thus, a technique 
for computing subsonic, oscillatory, blade-row air forces and moments 
would have several areas of application in the compressor field. More- 
over, the problem of tunnel-wall interference effects in subsonic wind- 
tunnel tests for oscillatory aerodynamic pressures constitutes a special 
case of the blade-row problem as treated herein, namely, the case in 
which stagger is zero and adjacent blades are 180° out of phase (this 
being the proper image system equivalent to the wall influence) . This 
problem has been previously treated by the National Advisory Committee 
for Aeronautics (refs. 1 and 2) by utilizing a form of aerodynamic inte- 
gral equation which constitutes a natural extension of the equation for 
the single wing oscillating in subsonic flow. It is believed that the 
computation method developed herein for the more general case including 
stagger and arbitrary interblade phase-lag angle may prove more conven- 
ient for application to the wall interference problem because of its 
relatively simple form of kernel function. It is even probable that 
the wall interference problem for the tunnel with slotted walls, whose 
image equivalent constitutes the special case of zero stagger and in- 
phase motion of all blades, may prove most easily calculable from the 
scheme derived herein. 

The present work constitutes a natural extension of previous inves- 
tigations (refs. 3 and 4) concerning the general problem of oscillatory 
blade-row behavior under flutter or vibration. These investigations 
established the blade -system mode shapes which characterize a blade row 
at flutter (or even in vibration distinct from flutter) and developed 
the computation techniques for oscillatory air forces in the incompres- 
sible flow regime. This report presents an entirely different approach 
to the problem of calculating aerodynamic forces and moments experienced 
by cascades or of approximating compressor blade rows in the subsonic 
two-dimensional regime. Considerations are limited to low incidence, 
and effects of blade thickness and camber are neglected. The mutual 
oscillatory aerodynamic interference of all blades is accounted for under 
the class of system mode shapes the existence of which is proven in 
reference 3* 

By use of the Fourier transform the kernel function of the aerody- 
namic integral equation is effectively reduced from an extremely cumber- 
some form involving semi -infinite -range integrals of doubly infinite 
series of Eankel functions to a simple, single, closed-form expression. 
Moreover, the nature of the kernel is such that close control over accu- 
racy is facilitated. Concurrently with the study of the formulation and 
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numerical solution of the aerodynamic integral equation considerable 
attention has been given to the occurrence of the so-called resonance 
phenomenon. This resonance phenomenon constitutes an infinite set of 
values of a parameter depending upon flow and configuration character- 
istics at which the aerodynamic kernel function becomes everywhere 
infinite. The phenomenon was described and treated in detail by Woolston 
and Runyan (ref. l) for the wind -tunnel -wall interference problem and, 
consequently, could be interpreted as a situation to be encountered by 
the cascade or blade row at zero stagger and under antiphase blade 
motion. The present report shows that the resonance condition may occur 
in staggered cascades with arbitrary interblade phase lag, the resonance 
parameter values depending upon Mach number, oscillatory frequency, 
interblade phase angle, gap, stagger, and acoustic velocity. 

The resonance phenomenon could be a matter of serious concern to 
compressor designers in view of the fact that it represents a set of 
conditions at which all self -induced aerodynamic forces, including the 
beneficial aerodynamic damping usually present below flutter speed, dis- 
appear. The blades act effectively as if they were in a vacuum. Thus, 
not only will certain critical blade oscillation frequencies be undamped, 
but disturbances caused by wakes or upstream disturbances giving rise to 
the proper combination of frequency and interblade phase lag may fall at 
the resonance condition with consequent disappearance of aerodynamic 
damping. This resonance -induced damping loss in an actual compressor 
configuration will, of course, be mitigated by the radial variation of 
cascade and flow parameters. These variations preclude the possibility 
of a major radial extent of the blade row becoming simultaneously resonant. 

Following the development of the aerodynamic integral" equation in 
the Fourier transform domain, the discussion of resonance, and the formu- 
lation of calculation techniques, the results of a specific set of com- 
putations are presented in order to provide confidence in the feasibility 
and accuracy of the method. The conditions for which the calculations 
are performed correspond to the conditions at which Woolston and Runyan 
(ref. l) performed computations for the wall interference problem by the 
more conventional but far more cumbersome Hankel function series type of 
kernel formulation. The agreement is excellent in all cases. The method 
described herein is thus believed to be "checked out" for its validity 
and feasibility and, thus, to constitute the most reasonable technique 
for computation of blade-row subsonic oscillatory air forces and moments. 

This investigation was carried out at New York University under the 
sponsorship and with the financial assistance of the MCA. The authors 
wish to express their gratitude to Messrs. Woolston and Runyan of the 
Langley Aeronautical Laboratory for their kind cooperation in making 
available the results of subsonic, oscillatory, single-wing calculations 
and for their encouraging attitude toward the present work. 
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SYMBOLS 


* * A ; 

( 0 ) A n^(i) A n 


a 

b 

C(k) 

c 


d* 


coefficients in pressure- jump expansion 

values of A n associated with constant and linearly- 
varying up wash, respectively 

dimensionless elastic-axis location 

semi chord 

Theodorsen flutter function 
chord 

free -stream sound speed 

slant distance between blades in Prandtl-Glauert space, 

# + (s*) 2 /p2 


F m = 3f m (a)/p 0 bU 2 


f m(a) 


Fourier transform of mth-blade reduced pressure jump 


Ho (2) < ) 

h 

j n ( ) 

k 

L 


L h ,L h 


L e> L e 


M 

^ea 


Hankel function, second kind, order zero 

translation, downward, of blade elastic axis 
Bessel function, order n 

reduced frequency, cbb/u 
lift amplitude 

lift amplitude due to translation; blade in cascade 
and isolated blade, respectively 

lift amplitude due to pitch; blade in cascade and iso- 
lated blade, respectively 

free -stream Mach number 

nose-up moment amplitude about elastic axis 
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w 


midchord moment amplitude due to translation; in cas- 
cade and isolated blade , respectively 




nose -up moment amplitude about midchord 


M e> M e ' 


midchord moment amplitude due to pitch; in cascade and 
isolated blade, respectively 


N 


magnitude used in replacing infinite integral limits 
by finite limits 


P 

M 

s 

s* 

t 

u 

u 

V 


pressure 

pressure jump, (p upper - P lower ) 

gap distance 
stagger distance 
time 

free-stream velocity (relative to blades) 
x-direction disturbance velocity 
y-direction disturbance velocity 

coefficients in the expansion of upwash for incompres- 
sible case 


x' = x/p 

j,y' 

CL ? T 

p = \fl - M 2 

7 

Tj = cn/c 0 p 

9 

A 


streamwise coordinates 


coordinates normal to stream 
Fourier transform variables 


angle utilized in angular transformation of chordwise 
coordinate 

pitch angle 

transformation coordinate defined in equation (Cl) 
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H = km/p 2 


£' = l/P 


r n ~ m 

r l ,t 2 


0 L A ' 

•41 4i 


<W 


= a 


aMs* 


Subscript: 


integers 0, ±1, ±2* . . . 

free -stream density 
interblade phase -lag angle 
poles of kernel function 

phase angle between translation and lift due to trans- 
lation; cascade and isolated blade , respectively 

phase angle between pitch and lift due to pitch 

phase angle between translation and resulting moment 
in cascade and isolated blade, respectively 

phase angle between pitch and moment due to pitch 
reduced pressure function defined in equation (6) 

frequency 

frequency of lowest resonance condition 

amplitude, under harmonic-time dependence 
real part 

imaginary part 

jump, (Upper - Lower) 
signum function 

corresponding to mth blade 



MCA TN 4136 


7 


FORMULATION OF PROBLEM 


A development of the linearized aerodynamics problem for the oscil- 
lating cascade in subsonic compressible flow is presented as follows. 

At the outset the system mode shape is introduced whose existence is 
established in references 3 and 4. This states, in effect, that at 
flutter the cascade or blade-row system oscillates with but a single 
characteristic blade mode, this mode occurring with equal amplitude in 
all blades and with a uniform phase shift between adjacent blades The 
ultimate value of this phase shift is, in the flutter of an actual 
rotor -mounted configuration, that value chosen from among a discrete 
set of admissible phase-shift values which minimizes flutter velocity. 

For a linear cascade and, practically, for a rotor -mounted blade row 
with many blades, it may be assumed that the phase angle is continuously 
variable and the minimization may be performed with respect to the 
resulting continuum. It should be noted at this point that the presence 
of resonance and its dependence upon phase angle in the compressible 
subsonic case might raise a question with respect to the validity of 
replacing the minimization problem for a discrete admissible set by a 
continuum-minimization problem. The fact that resonance will be auto- 
matically avoided by the system at flutter, however, indicates that the 
assumption of continuously variable phase angle, for purposes of flutter- 
velocity minimization in the many-bladed configuration, should remain 
valid in the subsonic case. In the event that the assumption of con- 
tinuously variable phase angle is rendered physically unjustifiable by 
a trend toward smaller blade numbers, then the minimization process 
must be performed with respect to the discrete set of admissible phase 
angles, equal to the number of blades in the blade row. 

The blade-row or cascade geometry is shown in figure 1, all vari- 
ables being physical or dimensional ones. Later scale lengths will be 
taken with respect to the semichord for convenience. The problem treated 
is two dimensional, initial blade angles are assumed to be zero or very 
small, and oscillatory disturbances are considered sufficiently small 
for linearized theory to apply. Under this theory the equation satisfied 
by velocity potential and acceleration potential, and hence also by pres- 
sure or either of the disturbance velocities, is 

(1 - + w - Si Wjrt - 1 w tt . 0 (1) 

0 c 0 

where W may be taken to mean any of the above-mentioned quantities. 
Writing the equation in terms of pressure and at the same time intro- 
ducing harmonic time dependence as 


iu>t 


p(x,y,t) = p(x,y)e 


( 2 ) 
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gives 


(1 - Maip^ + - 21 tt 2 P x + 7-5 5 = 0 


( 3 ) 


If the substitution 

x' = x/p 

y' = y 

is employed, equation ( 3 ) takes the form 


(4) 


■^x'x' + ^y'y' 


2iMa> ^ + of 

c Q P x ' c r 


P = 0 


( 5 ) 


The function \|r(x’,y’) is then introduced in accordance with the 
defining relation 


p(x',y') = iKx’ ,y' ) e in ^ x: / c oP 


In terms of i)f(x',y'), equation ( 5 ) in turn becomes 


tir + \lr + — ^ 

Vx' Vy' ~ 2 r 2 

c O P 


\|r = 0 


( 6 ) 


( 7 ) 


which is simply the familiar reduced wave equation. For convenience the 
parameter t] is introduced, where 


n = “/ c oP 


(8) 


Equation ( 7 ) then becomes 

Vx’ + Vy’ + ^ = 0 


( 9 ) 


At this point it is assumed that frequency co, hence also T) and later 
reduced frequency k = cbb/U, all have a small negative imaginary part. 


CD = CD-|_ + icD2 

k = k^ + ik2 l 

*1 = Tl]. + in 2 


( 10 ) 
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where 


0 < 


-a>2 

-k 2 


y« 1 


-T} 2 


The effect of this will he to replace the usual outgoing-wave require- 
ment at infinity by a requirement of boundedness at infinity. The 
device is a commonly employed one. Ultimately the imaginary part of 
a) will be forced to zero in a limiting operation. 


Singling out for inspection the mth blade, the following notation 
is introduced: 


y m ’ 

= r 

- ms = y - 

ms 

V 

- V 1 

ms* x 

> 

- ms* 

— IV 

P 

p 


( 11 ) 


Separation of variables shows that equation ( 9 ) has solutions of 

the form e^( ax '-y I n/ 'H -a ) which may be superimposed with respect to 
the separation parameter a. Then a solution to equation (9) with jump 
conditions in the proper form for the mth blade is 





( 12 ) 


where f m (a) is an undetermined function. It should be noted that \|r m 
as expressed in equation (12) satisfies the differential equation ( 9 ) 
(i.e., differentiation under the integral is permissible) as long as 
| y m ' | is held away from zero, provided the "proper" branch of the func- 
tion \j T] 2 - a?- is used and f m (a) behaves, at worst, algebraically 
at infinity. This proper branch is the one which behaves like -i|a| 
as | a 1 becomes large. Further, the slightly imaginary character of tj 
is needed to keep \|/ m bounded as |y m '| “ f° r | a l < | r li| • 

The jump |l|yj in the function \|r m from the upper to the lower 
side of the line y' = ms is given by 
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The fact that this jump must vanish for |x m '| > ^ will he taken care 

of subsequently. Moreover, the functions f m (a) will be shown ulti- 
mately to behave like Bess el functions of argument ab/p and hence to 
approach zero like l/\/|a| as |a| -»<». In view of this and chapter X, 
section 120, of Carslaw (ref. 5) > 'the integral in equation ( 13 ) exists. 

If the complete solution \|r to equation ( 9 ) is expressed simply 
as a sum over all blades of functions \|r m , then 

t = £ ( lJ+ ) 

m?=«oo 

Then it is evident that the jump [\|r J in \|r across the line y = ms 
is given by the term alone, all other values of \|r n being contin- 

uous across y = ms. Now from the phase-shift characteristic of the 
system mode shape under consideration, it follows that 

[%(*')] = - apjl = (15) 

where a is the angle. of phase shift between adjacent blades. Thus, 
from relation ( 6 ) there . follows 

[^(xOje^'AoP = e lmcr [ 4 r 0 (x m ')j e ;L ^ Mx:m, / Cof3 (l 6 a) 

or 

. / cuMs*\ 

im a ^ 

= [^(xm') e ' ' (l 6 b) 

In view of the fact the interblade phase -lag angle 0 will appear 
only in the form in which it enters expression (l 6 b), a new variable ft 
is defined as 

ft = 0 - 

At this point, it is assumed that a has a small imaginary part 0 2 , 

0 = 0- + i0 2 , with 0 2 = a> 2 Ms*/c 0 p 2 i n order to keep the quantity ft 
real. It should be recalled that ultimately o> 2 , hence also 0 2 , will 
be forced to zero. Henceforth, therefore, the quantity ft is to be 
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considered real . Just as the entire problem is periodic in a with 
period 2rt, so it will obviously be periodic in ft with period 2rt. 

Inserting equations (l6) into equation (13) and making use of equa- 
tion (17) gives 


/ OO 

^ f m (a)e iax «i' da 


(18) 


■ » Q 

Inverting this expression, considering JjijrJ e x and f m (a) as Fourier 

transforms of one another, and noting that jjr 0 (x')] = 0 for |x'| > — 
leads to the relation 


b/p 


*■<«> r [ to(% -)] 

^ Jy b/fi 


r ictx tn' 


x m ' — ’l*/ P 

= e imft f 0 (a) 


dx m ' 


(19) 


Thus, the first result of the system -mode -shape concept is that the 
Fourier transform f m (a) of the pressure-jump function over the mth 
blade is equal to that, of the pressure-jump function for the zeroth 
blade multiplied by e 


imft 


Utilizing equation (19) in equation (12) and combining all values 
of ^ , as in equation (l4), to form the complete expression for \|r, 


V = 


^imft sgn(y m ' ) 
2 


P“ 

I f 0 ( a )' 

vj «,oo 


i(ax m ' - |y m ' |Jri 2 -a 2 ) 

5 da 


(20) 


m=-°o 


Now, from the second Euler equation of motion, in linearized form, the 
y-derivative of pressure relates to the material derivative of upwash 
v as follows: 


v + + Uv x =~p. 


Po‘ y 


( 21 ) 


Under harmonic time dependence and in terms of ( x 
hecomes 


■ - F - y) 


this 


icov + - V Y ! 

p x 



■y 


t 


(22) 
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Equation (22) can be expressed in the form 


or 



(23a) 


(23b) 


In view of the fact that upwash must vanish infinitely far upstream, 
equations (23) may be integrated to give 


v(x’,y') . e-^'A r> 05 '/U/_ J_\p J V ,y) « . 

J -OO \ Pr^U/ ^ 


(24) 


Now expressions (6) and (20) may be utilized in equation (24) to get a 
relation between the upwash amplitude function v and the unknown 
Fourier transform of the zeroth-blade pressure-jump function as follows 



ms 

where | m ' = § ' — . Actually, the limiting operation lim^_ or 

lim _ or lim should be understood in equations (24) and (25 ) • 

® kg— » 0” 

Interchanging the order of integration in equation (25) , which is 
legal as long as |y m '| > 0, the expression for v becomes (utilizing 
the slightly negative imaginary component of cd) 


v(x' ,y ' ) 


= _J_ lim -tn|Bx'/U P* f / ) ~ a2 e l4Dx'/Up 

_ TT ... r\— I n v / , 


2p 0 U gv>— * 0 


Z imN il 
e e 


7 

(ax m '-|y m ' I \Zl 2 ^ 2 ) 


(° 




X 


m=-oo 


da 


(26) 
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v(x',y') = lim _ £ — — — 
032*-^ 0 2P Q U 


gl^x'/up P“ f Q (a) \/n 2 - a 2 


("ft) 


I 


imfl i(ax m , -|y m ' | , 


It is next verified that 


+ x',nsj = e incr v(x', 0 ) 

For, inserting jjx' + (ns*/p)] for x' and ns for y' in equa- 
tion ( 27 ) gives 


vfx' + ,n 

V 3 


rfx'/up e imM 2 ns*/UP 2 /*" f 0 ( a ) 


/ n 2 - a 2 e iax ’ da V e 1 * 10 e C °P e L\ P P 


, x r 

-j- 1 ns- ms | \J 1 


or, letting v = (m - n) 


v x- + M*,nsl = lim e ln<J 

V P ) 0“ 


fte U P r f fa (a)V n 2 - a 2 e iax> 




s W T) 2 -a 2 


is obtained as was to "be proven. This result is in accordance with the 
system-mode -shape concept and reduces the problem, by its periodic 
nature, from one in which boundary conditions must be inserted on all 
'blades to one in which they need be inserted on only one, say the zeroth 
blade. Therefore the upwash amplitude distribution v(x’) on the 
zeroth blade, which is assumed known for - ^ < x’ < is given by 


v<*',o). ii. r )E£Zl v 

a> 2 -> 0 ” 2 p Q U J -00 a + — ^ 

0 Up m=-«> 


e lax ' da, ( 51 ) 


The summation (bracketed expression) may be expressed, in closed 
form, as 


i s in ( s y rp - a?- 


,2 _ „2 U 


cos s wt)^ - ] - cos m 
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Utilizing equation (3l)> equation (30) becomes 


v(x',o) = v 0 (x') 


lim 

o>2-> 0“ 


ip e imMx7 c oP 

2^u 




(33) 


Once again the l/j/]"a| type behavior of F c (a) as |a|-»«> insures 
the existence of the integral expression for upwash. Now vq(x') is 
known only for | x'| < b/p. For |x'| > b/|3 the upwash amplitude func- 
tion vq(x') is unknown, but the zero-pressure-jump condition is known 
to hold. Thus, equation (32), considered as an integral equation for 
f 0 (a), constitutes a three -part mixed-boundary -value problem (ref. 6). 

This exceeds by one the number of parts amenable to Wiener -Hopf tech- 
nique. Had the blades been doubly infinite in the chordwise direction, 
then a straightforward use of Fourier transforms would allow solution 
of equation (33) for f 0 (ct) . For blades semi -infinite in extent, a 
Wiener -Hopf extension of Fourier transform technique would suffice. 

The actual problem, with finite blades, exceeds the capabilities of 
even Wiener -Hopf technique and hence demands some form of variational 
or numerical treatment. 


The approach to be outlined herein will thus constitute a numeri- 
cal technique for finding what amounts to the Fourier transform of the 
desired pressure function rather than the pressure function itself. 
This appears to be a somewhat novel feature of the method. As will be 
shown subsequently, it will be unnecessary to transform the solution 
function f 0 (a) (or the numerical approximation thereto) in order to 
get the pressure -distribution function. Rather, it will be shown that 
oscillatory lift and moment amplitudes are expressible directly in 
terms of f 0 (a) and its first derivative. 


The procedure is as follows: From equations (6) and (13), special- 

ized to the zeroth blade, 


p 0 (x'jl - [“♦ 0 (x')le teMx 'AoP = e taMx'/c 0 !3 f” f 0 ( a )e lc “' da (34) 
1 -00 


At this point, a new dimensionless system of variables is introduced. 
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x = ^l = 

b 

k = oob/U 
t = ob/|3 

Fh(t) = 


^ (Physical) 

D 


Pf 0 ( a ) 

77 ^ 


P f °(f 

PnW 2 


(55) 


Then equation (34) takes the form 

ijPo( x 0 ikM^x/p^ 


= e 


p o U 


/ 


~ / \ 1 TX , 

F Q (T)e dT 


while equation (33) becomes 


(36) 


!°<h, alm i fe i»Vs £ 

k^O" 2 


U 


f 

'J T~-a 


^ - t 2 PdCTje^Bin - t 2 


T + 


k \ 


COS 


Sp )k % 2 


W cLt 


b \| ^ 


- COS ( - 


TS 


* 


(37) 


Discussion of the first-order pole at t = -4L will be delayed until 

P 

the limiting process k]_-» 0“ is performed. 

It is interesting to note that, despite the presence of the term 


— in equation (37 ) > because of the manner in which the term 

. P 

appears, the entire integrand is single valued or possesses no branch 
points. This fact will not be used, however, and the original branch 
choice stated below equation (l2) will be followed. This is necessary 
because the integrand of equation (37) will subsequently be separated 
into two terms, each of which is itself multivalued or possesses branch 
points . 


LIFT, MOMENT, AND FORM OF TRANSFORM FUNCTION F q (t) 


Next the statement concerning the direct relation between lift, 
moment, and the function F g (t) is verified. From equation ( 3 6 ), it 
follows that 
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( T ) = _L P 1 [goMi e -ikM 2 x/p 2 e -irx 

° 2n J-l Pq u 2 


dx 


(38) 


Thus, the dimensionless oscillatory lift amplitude is given by 


PnU 2 b 


= -2rtFJ- 


kM c 


°l p 2 


(39) 


considered positive upward, as in figure 1. Moreover, nose-up moment 
amplitude about midchord easily follows also from equation ( 38 ) . 


Mq 

P 0 U 2 b 2 


2*IFo 



(^0) 


Thus a solution for the Fourier transform F 0 (t) constitutes, in effect, 
a complete solution to the problem. 

The aerodynamics problem now reduces to the determination of a 
function F 0 (r) satisfying equation (37) with v Q (x) given in (-1,1) 
and such that |j) 0 (x)] vanishes for jxj >1. Referring to equation ( 36 ) 
for [j? 0 (x)j and to the treatment of Weber's discontinuous integrals in 
reference 7 it may "be seen that an expansion of F q (t) in a series of 
Bessel functions Jn(r) would automatically satisfy the zero -pres sure- 
jump condition off the blade proper (i.e., for | x| > l) Moreover, in 
view of the expansion formula (ref. 8) for the Bessel functions of a 
sum-type argument it is evident that an expansion of F 0 (t) in a series 
of Bessel functions J n (r + c), where c is a parameter independent of 
t, would serve the same purpose. Thus, by assuming F 0 (t) to be given 
by an expansion in Bessel functions Jn(T) the off -blade-pressure con- 
tinuity condition is automatically satisfied. All that remains is to 
determine the coefficients in the expansion such that equation (37) is 
satisfied to the required degree of approximation on the wing (i.e., 
for | x| < l) . 

In order to facilitate the solution to the problem and, at the same 
time, put F q (t) in a form such that behavior of the integrand in equa- 
tion (37) may "be examined, a trigonometric transformation together with 
a corresponding pressure expansion, common in both steady and nonsteady 
two-dimensional wing theory, is introduced as follows: 
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X = —COS 7 


P Q (x) 

= A 0 cot 

P 0 U 2 


00 


An sin ny 


(M) 


Expression (Al) for pressure jump contains the proper leading-edge 
singularity and vanishing trailing -edge value of pressure jump, in 
accordance with one form of the Kutta condition. Substituting equa- 
tion (kl) into equation (38) leads to an expansion form for F q (t): 



Wow, as could be verified by direct integration of equation (kl), the 
dimensionless lift amplitude is given simply by referring to equa- 

2 

tion (39) and letting t = — ^ . Noting that the zeroth-order Bessel 

P 2 

function is the only contributor, there results 


L 



(43) 


Again, referring back to equation (kO) and utilizing the Bessel function 
recurrence properties, the moment is found to be given by the following 
expression: 


Mp 

P 0 U 2 b 2 



+ 



Ji'(o) 



(kk) 


Once more, this is easily verified by direct weighted integration of 
equation (kl) . 

The next important point to notice is that, as a result of the 
assumed expansion equation (kl) of the pressure- jump function, F q (t) 
takes a form, involving Bessel functions, which insures the vanishing 
of off -blade pressure jump. Thus there remains only to calculate A q, 
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A^_, and Ag. Unfortunately, however, as is usually the case in such 

problems, the higher index coefficients influence these required first 
three coefficients . 


Singularities and Resonance Condition 

It should be noted that the integral of equation (37) may possess, 
in addition to the first-order pole at t = ■=)§ (which approaches the 

3 2 

real axis as kg 0") , poles at the zeros of the cosine -difference term 
in the denominator. The form of this term makes it evident that these 
poles, if they lie on the real axis as kg -» 0”, must lie at values of 

k x M 

t so that | t| < It will now be shown that the arrangement of these 

1 1 p 2 

poles relates to the criterion for resonance. This is best accomplished 
by expressing the cosine -difference term in a different form 



-2 sin! 




sR jk 2 M 2 

pv 3^ 



( 45 ) 


Letting 

M 2 = k 2 M 2 /(3^ 

(a*) 2 = s 2 + hh! 

P 2 


(46) 


(Note that d* is the slant distance between blades after a Prandtl- 
Glauert transformation of the physical plane.) The first sine factor 
vanishes whenever 


IP UF 


2b 


O TS* 

= itn + 

2 2b 


for any integer value of n, positive, negative, 
sine factor vanishes whenever 


TS' 


\||i 2 - t 2 = jtm + — 

2b v 2 2b 


(47a) 

or zero} the second 

(47b) 
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for integer values of m. Calling the first set of poles T-^ n and the 
second set Tq 111 , solution of equations (47) gives 


T 


n 


1 


bs* 

p 2 (d*) 2 


(ft 


2«n) + 


(d*) 2 


V 


P 2 (d*)‘ 


— (ft - 2jtn) 2 

P 2 


(48a) 



— s *— (0 
p 2 (d*) 2 


+ 


2jon) 



(ft + 2«m) 2 

P 


(48b) 


where the poles T-j_ n , when complex, lie in the lower half plane and the 
poles T 2 m , when complex, lie in the upper half plane. It is of interest 
to note that, as |n| or |m| becomes large, 



-2rtb /ns* 

p(d*) 2 V P 


1 n s 


(49a) 


T m _> _2«b pSl + i| m | a \ (4 9 b) 

2 (S(a*) 2 l p I 

In other words, for large values of |n| and |m| the poles approach 
positions similar to the cascade itself together with its reflection in 
the y-axis after a Prandtl-Glauert type of transformation. Unfortunately, 
the possible utilization of residue calculus in the evaluation of the 
infinite integrals such as appear in equation (37) is defeated by the 
behavior at infinity of the function F 0 (r)ei' rx for | x| < 1. Since 
F 0 (t) is known to be expressible in an expansion of Bessel functions, 
it, therefore, acts like cos r/\J~T for large values of t, and 
F 0 (-r)e lTX fails to converge on either an upper or lower semicircular 

contour (for | x| < l) as the contour radius approaches infinity. The 
impossibility of evaluating the infinite integrals via the residue cal- 
culus reduces the behavior of the complex poles T-^ n and r^°- to a 

matter of academic interest only. The numerical solution of the inte- 
gral equation must be performed with integrations over the real line 
as they originally occur. Thus the remainder of the analysis will be 
confined to real -line integrations . 


The existence and location of real poles T-j_ n and t ^ under cer- 
tain conditions may be related to the resonance phenomenon in the fol- 
lowing manner: For sufficiently small values of i_i 2 (d*) 2 p 2 /b 2 = 

k 2 M 2 (d*) 2 /b 2 p 2 there will be a range of values of Q. (recall that only 
values of ft in the interval (0,2rt) need be given consideration because 
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of the periodicity of the entire problem in ft) in which ft is bounded 
away from all multiples of 2k (including zero) far enough to prevent 
the radicands in equation (48) from becoming positive. For such com- 
binations of ft, k, M, and d*/b no real roots are possible. This 
situation corresponds to the region labelled A in figure 2. As noted 
in the figure, no real poles of the "cascade term" exist for this zone 
which will be shown to be below the "lowest resonance point." Now as 
either of the two boundary lines of zone A are crossed one real root 
T^ 11 appears, and simultaneously a real root T2~ n appears. Moreover, 

they first appear at precisely the same point on the real r-axis, indi- 
cating the presence of a second-order pole. This second-order pole 
constitutes a nonintegrable singularity and hence suggests a resonance 
condition. As shown by figure 2, the condition occurs at 

ft = kMd*/pb (50a) 

or 

ft = 2k - (kMd*/pb) (50b) 

these being the boundaries of zone A. As zones B and C are entered, the 
real poles persist but separate into two first-order poles. They thus 
constitute a pair of Cauchy type singularities which are integrable in 
the sense of Cauchy principal values. Next, the boundaries of zone D 
are crossed, and two more poles appear on the real axis. Once again, 
they first appear simultaneously at the same 
order pole indicating resonance conditions, 
zones C and D is characterized by 

n . Ml 

|3b 

while that between zones B and D corresponds 

ft = 2rt - Mfl* 

Pb 


point and form a second- 
The boundary between 


(51) 
to 

(52) 


Entering into zone D, four poles exist on the real axis, but they are 
distinct and thus integrable in the Cauchy principal value sense. Pro- 


ceeding, in this manner, to successively higher values of more 

and more new poles appear on the real t-axis and, as they first appear, 
they do so in pairs, the pair first occurring at a common point as a 
zone boundary is crossed. The general equation for the zone boundaries 
is given by 


kMd*/pb = ±(2jtv - ft) 


(53a) 
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or 


0id*/c o p = ±(2jcv - 0 .) ( 531 ) 

with v an integer, positive, negative, or zero. 

Now condition (53) is precisely the set of resonance conditions for 
a staggered cascade at arbitrary interblade phase angle. For, when the 
present problem is formulated in terms of the Hankel function series, 
as given in appendix A, which is then converted via the Poisson summation 
formula (ref. 9) to an exponential series in much the same way as is 
performed in reference 1 for the case of tunnel-wall interference, then 
the resonance conditions appear as relations (53) • As is to be expected, 
these resonance conditions reduce to those of Runyan, Woolston, and 
Rainey (ref. 2) for the wall -interference problem upon specialization 
of d* to s (s* = 0) and fi (Q = a when s* = 0) to it. Since 
resonance, as herein defined, corresponds to a divergent kernel in 
equation (Al), for all values of x 0 , it follows that the pressure jump 
itself must vanish at the zone boundaries of figure 2. These boundaries 
are therefore resonance loci for the cascade problem. For the present, 
the method of computation of ' aerodynamic oscillatory lift and moment 
coefficients is confined to zone A, the region lying below the lowest 
resonance condition. 


Computation Methods 


The fundamental integral equation (37) is now treated for subreso- 
nant conditions (zone A). Referring to equations (31) and (32) it 
becomes evident that, in the case of a single wing oscillating in sub- 
sonic flow, the cascade term 


i sin s \/ T) 2 - a 2 j 

cos s \j -a 2 - cos - 2^ 

would be replaced simply by 1 (unity) . For purposes of computation 
single -wing term is added and subtracted in equation (37) "to give 


( 5 ^) 


, the 


V Q (x) 0 1 kM 2 y /ftS r 00 i-rjr 

0 Single -wing contribution + lim - 0 / p / dT - V 

/ k 2 M 2 2 

- 7 

i 


lA 2 _2 

U " ““ ~~ ' kg-KT 2 J t=-« ' ' T 


COS ^1, 

[ *\l 

- t 2 - cosiin - in') 

K* \ r> J 


(55) 


The reasons for handling the problem in this way will become evident as 
the development progresses. The first point to notice is that the brack- 
eted term in equation (55) converges exponentially to zero, thereby 
effectively reducing the li m its of integration from (-oo^co) to quite small 
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finite values for practical computational purposes. In fact if the infi- 
nite lower and upper limits are replaced by -A and B, respectively, then 
the error committed by neglecting the "tails" of the infinite integral 
may be shown to be bounded by e, where 


e 


0.4b 


1 J3s 
R b 


ei- 



2 


BP 


psR 

b 



V 


R 2 - 


k 2 M 2 

r4 



(56) 


and where R is the smaller of A and B. For reasonable values of 

M, — , and k, extremely high accuracy is obtainable with rather small 
b 

limit magnitude R. For example, with 

s* _ s _ o 
b b 


M = 0.5 
k = 1 


it is easily possible to have a subresonant range of fl, and a value of 
R = 8 gives an error of less than 1 x 10“°. 

Thus, expression (55) is replaced by a numerically equivalent expres- 
sion in which the finite limits -A and B replace the infinite limits. 

In particular, it proves convenient to take 


A = 



B = N 



(57) 


where N is such that the smaller magnitude N - 


_k_ 

„2 


must satisfy the 


error requirement (eq. (56)) on R for any prescribed accuracy e. As 
k/p2 will seldom in practice exceed 1.5 for M = 0.5 , it is seen that 
N = 9-5 will be consistent with R = 8.0 in the above example. With 
the limits now centered about the singularity t = — , the singularity 

p2 


is easily removed by the artifice of adding and subtracting the value 
of the remaining portion of the integrand, evaluated at t = Thus, 
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v D (x) .. . R ikM 2 x/p 2 ~ dr 

— = Single-wing contribution + lun £ e /H H EL 

U k 2 H k l + - 




G(t) = e iTX F q (t) 



i Binf& - 

b Vl^T 


s ft i /k%^ 1i 


b V p 4 


T - COS 


( n - I r) 


At this point the limiting operation k^-* 0“ is introduced. This has 
no untoward effect on the bracketed term in view of the continuous behav- 

[ G(T) ' “(#)] ] ,2 , 
ior of the quotient < J= — - ^ \ as t approaches -k/p^ (i.e., 


ior of the quotient 1 l 1 

as t approaches -k/p^ (i.e., 

1 ( T + i&) I 

in this term k is merely replaced by kq) . The expression in brackets 
g/ 1^§\ in equation (58), however, contributes a term outside the inte- 

• VP 2 / 

gral as the limiting operation is introduced 




k 2^°" -N- 


k l IT + 


L dT = 

.2 


4-h a l±i 


K-i / 

-N - -i T + -1 


dr (60) 


In equation (60), because of the symmetry of the interval of integratioi 

-k 

about the point t = and of the antisymmetry of the function — - — 

3 2 , 1 

f 

in the same range, the integral vanishes and only the external term 
remains. This explains the reason for symmetrizing the range of inte- 
gration in equation (57) • 

Performing these operations and replacing kq by k, assumed real 
henceforth, gives the equation 
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v ( x ) in i kM 2 x/p 2 , k > fi iiM 2 X^ 2 / 

L i Single-wing contribution + ^ pe \~) + 2 6 / 


i 5 




dr 


(6l) 


Now 


G 


^ = * lkx/p 


;ik\ 

0 


sinh 


sk 


cosh ^ - cos (a + TS2?\ 


- 1 


lip" 


u 


Hence ^ equation ( 6 l) becomes 


VpU) 

u 


/-tt-lkx 1 

sinh £ 

... 

P e ikM 2 x/P 2 

f-t 

[ Q(T) - G (il)l 

w 

ks / 

cosh — - cos (SI 

\ 

* >“l\ 
3 2 *j 

s J 

-n-Jl 

3 2 

( T + ?) 


(62) 


dT (63) 


for |xj < 1 where it is no longer necessary to indicate the Cauchy 
principal value since the integrand is now continuous at t = -k/p^ 
and where 



i[A n + ^Vi(- k ) + \ 



(64) 


and where G(t) is given by equation ( 5 9)* 

The form of equation ( 63 ) is the one recommended for computational 
purposes . It is to be emphasized that it holds only for the subresonant 
condition, since when poles t-J 1 and t g m lie on the real axis as 

kg -» 0 ” (superresonant conditions), then other terms are carried out- 
side the integral in the performance of the limiting operation indicated 
by equation ( 60 ) . 

Equation ( 62 ) is to be solved for the function F d (t) or, more 
precisely, for the first three coefficients Aq, A-j_, Ag, in the expan- 
sion (eq. (42)) of F g (t). Several methods come to mind, but the pres- 
ence of the single-wing contribution restricts the technique to that in 
which this contribution is expressible in the required form. In any 
case, under the assumed expansion (eq. (4l)) of the pressure -jump dis- 
tribution, methods must be used in which the contribution of the single 
wing to the upwash distribution is expressible in terms of the coeffi- 
cients A n of the pressure series. Two means of solving the resulting 
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equation are the pointwise or collocation method and the termwise or 
Galerkin method. In the former, the desired relation (eq. (62)) is 
satisfied at a finite number n of discrete points x in the interval 
(-1,1) and the same number n of coefficients is retained in the pres- 
sure series (eq. ( 4 l)). The resulting equations are satisfied for the 
first n coefficients A m . While only the lowest three values of A 

are required, nevertheless as more values of A are included in the 
collocation calculation, the accuracy with which these first three coef- 
ficients are calculated increases . This is the method recommended for 
the subsonic compressible case (0 < M < l) in view of the existence of 
numerical data obtained by Woolston and Runyan (ref. l) for the contribu- 
tion at each of three points x (-0.5, 0, and 0.5) to the upwash from 
the first three terms in a pressure series equivalent to expression ( 4 l) . 
In the incompressible case, a Galerkin technique is recommended, in which 
the error in equation (63) Is orthogonalized with respect to a finite 
number of the complete set of functions 1, cos 7, cos 2 7, cos 3 7 > 

. . . after replacement of x by (-cos 7) throughout equation (63). 

This procedure is outlined in detail in appendix B. 


COMPUTATION METHODS IN SUPERRESONANT CONDITION 


Returning to formulation ( 55 ) of the aerodynamic integral equation, 

considerations are extended to the case where zones to the right of the 

subresonant zone A (fig. 2 ) are entered. Hence real-axis singularities 

appear because of the vanishing of the denominator expression ( 45 ) in 

-k - 

addition to the obvious singularity at t = — ^ in the limit as ko -» 0 . 

0 2 

In each of zones B and C, there will be a total of three singularities, 
and in zone D there will be a total of five, and so on. The computations 
must, of course, be performed with parameter combinations lying within 
the interior of a zone, that is, bounded away from the zone boundaries 
or resonance lines. Hence, in the limit as k 2 -» 0 “, the real poles 

will be separated by nonvanishing distances. Consider the general case 
of n real poles t v (v = 1, 2, . . . n) (in the limit as ^ -* 0”), 

lr- 

one of these being the pole at r = — — and the others being the terms 

P 2 

T^ n and t ^ 1 of equations ( 48 ). Rewriting equation ( 55 ), f° r conven- 
ience^ in the form 

= Single-wing contribution + lim F h(t) (It (65) 

^ — > 0 ^ t= — N 
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where a value of N is chosen which is sufficiently large to give the 
required accuracy, centered intervals I v are introduced about each 
pole t v as shown below: 



The distinctness of the poles r v insures that nonoverlapping intervals 
I v may be constructed, although this is not crucial to the argument. 
Outside of the intervals I v but within the integration range (-N,N) , 
the cascade contribution (the integral contribution) of equation (6^) 
is calculated by numerical integration in straightforward fashion. 
Within an interval I v centered about the pole t v , the following pro- 
cedure is followed: 

The integrand H(t) may be written 

H(t) = H(t)(t - t v ) ^ j — 'j 

■ j(t ' t '' ) (tvt ;) 


= J(t>t v ) - J(t v ,t v ) + J(t v ,t v ) 


(66) 


T - T. 


T - T, 


where J(t,t v ) = H(t)(t - t v ) is regular at t = t v , since t 


is a 


simple pole of H(t). Then, in the limit, as k 2 -» 0~, the contribution 
of the term 


J(t,t v ) - J(t v ,t v ) 
t - t v 

to the interval I v of the integration is simply 

J(t,t v ) - J(t v ,t v ) 


(67) 


I, 


T - T, 


dr 


( 68 ) 


with real values of t v being understood. As k 2 -» 0“, the term 


T V ?T V ) 


of equation ( 65 ) contributes nothing internally to the inte- 

r - T V 

gral in view of the symmetry of the interval and the antisymmetry of 

the function — 1 , but a term 

t - r v 
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iirt j( Tyj Ty) (69) 

appears externally (ref. 10) (outside the integral), the plus sign 
corresponding to a pole t v which approaches the real axis from above, 
and the minus sign corresponding to approach from below. The pole at 

t = — is a special example of the first case, since as ko -» 0 “, 

P 2 

t = — approaches the real axis from above thereby contributing a term 
13 

(TO) 

as was shown in equation (60). 

In this fashion, the superresonant cases may be computed by a 
routine though somewhat more cumbersome procedure. 

It is to be noted that for the pole at t = -k/p^ the operations 
indicated in equation (66) are actually unnecessary in view of the fact 

that this singularity is already in the form Nevertheless, 

t - T V 

for the sake of uniformity, it does no harm to include this pole in the 
general considerations of equation (66) . 

It will be necessary to know the value of J(t v >t v ) in the limit 
as k2 -> 0“ or as the values of t v become real. Referring to equa- 
tion ( 55 ) for H(t) and to expressions (45) through (49) for the values 
of t v , the following steps may be taken: 



Hence, 

J(t v ,t v ) « lim H(t)(t - t v ) 
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In the case where t v is one of the r-j_ n terms of equation (48) with 
approach to the real axis from below, as k 2 — > 0", there results 


J(t v 


v) “ I 





(73) 


with x y 2 = JL ^2jtn - ft + - j and with t v 

tion (48a) . Similarly when t v is one of the terms 
tion (48b), approaching the real axis from above, 

F Q (T v )(p2 - t v 2) / 

( T v + \ T v s P “ 

with \j \x 2 _ Tv 2 = JL ^2jtm + ft - — j and with t v 

equation (48b) . 


J(' 


v> 1 v 


) = be 


IX 


given by equa- 
T 2 m of equa- 



given by 


RESULTS OF SAMPLE CALCULATION FOR ZERO STAGGER, 
ANTIPHASE CONDITION 


In order to establish confidence in the Fourier transform solution 
technique described herein, a sample set of calculations was carried out 
completely for a particular condition at which previous results, both 
experimental and theoretical, exist, the previous theory being that 
associated with the kernel function in the form of an infinite integral 
of an infinite series of Hankel functions (appendix A) . Runyan, Woolston, 
and Rainey (ref. 2) present experimental wind -tunnel-wall interference 
effects on oscillatory lift and moment in pure pitching motion (about 
midchord) for a series of subsonic Mach numbers with a ratio of tunnel 
height to wing chord of 3.80. They provide theoretical results, obtained 
from a three -point collocation treatment of the aerodynamic integral 
equation with inf inite -integral— inf inite -series kernel, for these same 
conditions as well as for pure translatory motion and for other tunnel 
heights. The tunnel-wall interference problem is, of course, completely 
equivalent to the blade-row problem under the restriction of zero stagger 
(s* = 0) and antiphase motion (c = 180°) . 
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The present Fourier transform technique, with its simpler, closed- 
form kernel, was utilized in a three-point collocation solution for the 
particular case of a blade row at zero stagger under antiphase motion, 
with a gap-chord ratio of 3-80 and for the particular value of Mach 
number M = 0.5* For this value of Mach number and for the range of 
values of reduced frequency (0 < k % 0.7) considered in the present 
example, three-point collocation gives satisfactory accuracy as has 
been proven by comparison with results obtained from five and seven- 
point collocation schemes for the isolated wing. The present computa- 
tion procedure requires the availability of single-wing or isolated- 
blade contribution to upwash resulting from each term A n in the 
pres sure -jump expansion as well as the contribution from the remaining 
blades in the blade row (or the images in the wall -interference inter- 
pretation) . This single-wing contribution was kindly provided by the 
Langley Aeronautical Laboratory of the MCA in the form of tables of 
integrals of the nonsingular portion of the single -wing kernel function 
with weighting functions corresponding to terms in a pressure-series 
expansion equivalent to equation (4l) herein. The contributions of the 
singular (logarithmic and first-order pole) terms were calculated and 
added to the MCA tabular values to provide the data of table I pre- 
sented herein. This is the tabulation of total single-wing contribution 
to upwash at three stations ( quarter -chord, midchord, and three -quarter - 
chord) due to each of the first three terms in the pressure-jump expan- 
sion (eq. (4l)). As noted in the table, the values are for M = 0.5 
and reduced frequency values k of 0.02, 0.1, 0.2, 0.3> 0.4, 0.5, 0.6, 
and 0.7- In the present case, reference to figure 2 and a brief calcu- 


lation show that the lowest critical value of k is k 

cnt 


0.715973. 


Table II presents the results of the computation scheme outlined in 
appendix C for all the remaining blades of the blade row, that is, for 
the terms exclusive of the single-wing contribution on the right-hand 
side of equation ( 63 ) . The corresponding values from tables I and II 
are simply added to obtain the upwash contribution of each pressure term 
for the complete blade row or the complete wing -plus -image system in the 
wall interference interpretation . 


Tables III and IV give the results for the single wing and the com- 
plete blade row, respectively, of the collocation process where the coef- 


ficients . A„ 

(0) n 


correspond to a constant unit relative upwash 



and the coefficients 


(l) A n 


correspond to a linear upwash variation 


( — = x ) inserted into the left-hand side of equation ( 63 ) . The expres- 

\ u ) 


sions below tables III and IV give the dimensionless lift and moment 
(complex) amplitudes as linear combinations of these A coefficients 
for arbitrary elastic -axis location and for arbitrary combinations of 
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pitching and translatory motion. The manner in which the systems (q)^ 11 
and are utilized in computing air forces and moments leads to 

use of the term "influence coefficients" in describing them. 


Figures 3 to 10 illustrate the comparison of lift and moment coef- 
ficients derived from the present theory with both the theoretical and 
experimental results of Runyan, Woolston, and Rainey (ref. 2). These 
figures relate to either pure pitching motion about midchord or pure 
translatory motion. It should be emphasized, however, that the data of 
tables III and IV permit simple calculation of lift and moment for any 
combination of translation and pitching about any axis. Figure 3 gives 
the results, in the form of lift amplitude ratio for the blades in 
cascade or the wing between walls to that of the isolated wing, for 
pure pitching motion about midchord. The present theory denoted NYU 
theory is seen to follow closely the theory of reference 2 denoted 
MCA theory, whereas both theories depart noticeably from the experi- 
mental results. The major portion of this discrepancy between theory 
and experiment probably lies in the pressure-integration technique 
utilized in reference 2 to obtain lift (and moment) from chordwise 
pressure -jump measurements. Figure 4 gives theoretical and experimental 
phase angles for lift in pure pitching about midchord, and once again 
the NYU and MCA theories follow each other quite closely. The single- 
wing curve, as shown for purposes of comparison in figure 4, can be 
obtained from the available single-wing tables or from a simple calcu- 
lation using table III. Figures 5 and 6 give amplitude ratios and 
phase angles for moment about midchord in pure pitching about midchord. 

The two theories give very close agreement, and the experimental ampli- 
tudes lie at a nearly constant distance below the theoretical amplitudes, 
again probably because of the pressure-integration method of obtaining 
experimental moments . The MCA experimental values and NYU theoretical 
values of moment phase angle are seen to correspond closely in figure 6 
up to just below resonance. As explained in reference 2, the pressure 
integration technique used therein would be expected to exhibit less 
severe effects on phase angles than on amplitudes. 

Figures 7 to 10 relate to pure translatory wing or blade motion for 
which MCA theory but no experimental data are available. These figures 
thus provide simply a comparison of numerical results from the two theo- 
ries. The agreement is seen in all cases to be excellent, to such a 
degree, in fact, that the curves are in some regions indistinguishable. 

If a generalization were to be attempted regarding the comparison 
of NYU and MCA theoretical results, it might be to the effect that the 
NYU results emphasize slightly the "saddleback" shape characteristic of 
the amplitude -ratio curves and exhibit a sharper decrease in phase angle 
as the resonant frequency is approached. On the whole, however, the two 
theories give extremely close results and, thus, each affords an increased 
measure of confidence in the other. 
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REMARKS CONCERNING APPLICATION OF METHOD TO 
BLADE -ROW PROBLEMS 


The extremely good agreement between MCA calculations based upon 
the form of kernel function having an infinite integral and infinite 
Hahkel function series and the NYU results based upon the simple closed- 
form kernel function stemming from the use of Fourier transform serves 
to justify the steps taken in the latter theory and to give this method 
a "clean bill of health" so that it may be used with confidence in the 
more cumbersome blade -row problems involving nonzero stagger and inter - 
blade phase-lag angles other than l80°. Comparison of equation (63) of 
the present report, which includes the general case of arbitrary stagger 
and arbitrary interblade phase-lag, with expressions in reference 2 for 
the kernel function in the more restrictive case of zero stagger and 
antiphase motion serves to illustrate the marked simplicity achieved by 
the transition to Fourier transform and the subsequent summing of the 
blade-row series and integration of the semi-infinite integral necessary 
to convert from pressure to upwash. Furthermore, the exponential decay 
of the integrand of equation (63) permits setting the finite integration 
limits to achieve any accuracy desired so long as a quadrature scheme 
of compatible accuracy is utilized. The use of a three-point colloca- 
tion technique probably constitutes the weakest link in the accuracy of 
the present calculations; however, five- or seven-point schemes would 
require the availability of five- or seven-point tables for the single, 
isolated wing. 


For application to turbine or compressor blade-row flutter or 
response problems the aerodynamic influence coefficients ( 0 ) An an< ^ 


(l) A n 


must be computed for an appropriate range of Mach numbers and 


reduced frequencies and for a range of values of D between 0 and 2n 
(this range corresponds to a shifted range of length 2it for the inter- 
blade phase angle 0 at any particular combination of the parameters 
o>, M, s*, and c Q ) . Critical flutter conditions then correspond to 
the minimization of flutter velocity with respect to 0 as shown in 
reference 3* Actually, as explained in that reference 0 must vary, 
not continuously in the range 0 to 2 jt but rather among a discrete set 
of values depending upon the number of blades in the blade row. How- 
ever, for large blade numbers, the assumption of continuously varying 
values of 0 is shown in reference 3 to be justified in facilitating 
the minimization procedure. 


The emphasis placed herein on subresonant conditions (fig. 2, 
zone A) is justified for configurations and operating conditions under 

kMd* 

which expected values of the parameter — do not exceed a value of 
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unity during the occurrence of flutter. For values of this parameter 
less than unity, critical conditions will lie in a vertical strip adja- 
cent to the ft -axis in figure 2, the strip width being less than one- 
third of the altitude of the zone A triangle. While this in itself 
does not insure that flutter cannot occur in a superresonant condition 
(zones B or C), all indications are that it does not. Calculations 
performed by Lane for the incompressible case M = 0 have always indi- 
cated a distinct boundedness of critical flutter conditions away from 
interblade phase angles a, of 0 and 2n. In the compressible case, the 
vanishing of self -excitation air forces along the resonance lines aa 
and bb (fig. 2) precludes the possibility of. flutter, in the immediate 
vicinity of these lines, of a blade row having any internal damping. 
Hence all indications are that blade-row flutter would occur with values 
of Q (note Q rather than a in the compressible case.) between 
(l/4) Jt and (7/4 )t: and with critical conditions falling in the subreso- 
nant zone A. When values of the parameter in excess of unity are 
expected, then calculations in the superresonant zones B, C, and D are 
necessary. The necessity for performing computations in still higher 
zones (E, F, etc.) appears most unlikely in any case. 


CONCLUDING REMARKS 


The objective of the work reported herein was to develop a feasible 
computation technique for the aerodynamic oscillatory lifts and moments 
experienced by cascades or compressor blade-rows in subsonic flow as well 
as to examine analytically the existence of the resonance phenomenon and 
the location or set of conditions at which the kernel function of the 
aerodynamic integral equation blows up. The objective was achieved by 
resort to a new technique in which the aerodynamics problem is formulated 
in terms of the Fourier transform of the blade pressure distributions and 
numerical solutions are found for the transforms rather than for the 
pressure distributions proper. It is shown that oscillatory lifts and 
moments are easily derivable from the transforms without the necessity 
of inverting. Moreover, the technique permits the expression of the 
kernel function in a simple closed form, whereas conventional methods 
lead to an extremely cumbersome kernel consisting of a semi -infinite 
range integral of a doubly infinite series of Hankel functions of the 
type encountered in single-wing subsonic oscillatory aerodynamic 
considerations . 

A plot is presented giving the location of resonance conditions, 
in general, and calculation results are reported for a configuration 
corresponding to the tunnel-wall interference problem which constitutes 
a special case of the cascade problem. The conditions for which calcu- 
lations were performed correspond to the conditions at which the National 
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Advisory Committee for Aeronautics had previously executed wall inter- 
ference computations following the conventional integral-equation formu- 
lation. The agreement achieved is excellent. The technique developed 
herein is thus believed to represent a valid and feasible method for 
obtaining oscillatory cascade or blade-row aerodynamic forces in the 
subsonic regime. 


New York University, 

New York, N. Y., August 2, 1956. 
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APPENDIX A 

EXPRESSION OF AERODYNAMIC INTEGRAL EQUATION 
IN CONVENTIONAL FORM 


For the sake of completeness, the integral equation for the 
pressure-jump amplitude distribution is presented below in a more con- 
ventional form, that is, the form involving a kernel expressed as an 
infinite series of Bankel functions: 

VpU) = 4 y e imfl r 1 e “ik(x— X q) [?q(^o)] r x - x o e ikr)/p 2 d£_ ^(2) 

D Lj Jx^-l PqII 2 j !)=-■» dy 2 

m=-oo 

where y, rj, x, and x D are all scaled with respect to the semichord 
b, that is, all are dimensionless. The above integral -equation formu- 
lation results from a straightforward application of Green's function 
for the reduced wave equation to the cascade or blade-row problem tinder 
the type of blade -to-blade periodicity predicted by the system mode 
theory. Once again, the interblade phase-lag angle a appears only in 

( cdMs - ^* \ / \ 

a ] defined as ft in equation ( 17 ) • 

. °oP 2 / 




di) dx 0 


(Al) 
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APPENDIX B 

APPLICATION OF AERODYNAMIC INTEGRAL EQUATION 
TO INCOMPRESSIBLE FLOW 


The aerodynamic integral equation, equation ( 65 ), becomes, in the 
case of incompressible flow (c 0 -» °°) , 


^ 2.-.— = Single-wing contribution + — F (-k)e 
U 2 ° 


sinhl 


cosh ~ - cos ^cr + 


- 1 


1 r N ' k Oq(t) - gp(-k) 
^ w -N-k ( t + k) 


cLt (B] 


noting the fact that, as c 0 00 , ft -» a. In equation (Bl) G 0 (t) 

signifies the limit of G( t) as c 0 -» <» and is given by 


G 0 (t) 


-1 t e 


1TX 


■F 0 (t) 


sinh 

ST 


b 


cosh/— N ] - cos 1 
\ b 


:-?) 


- 1 


(B2) 


Equation (Bl) holds for all values of cr, k, s, and s* since the 
possibility of resonance disappears when c Q -» 00 , This is evident from 

examination of figure 2 where, with zero Mach number, the entire problem 
degenerates to the ft-axis and only the points ft = 0 and ft = 2it (or 
a = 0 and a = 2jt) appear as resonance conditions. Upon closer examina- 
tion, however, it becomes clear that even these two points do not give 
resonance. This is because, with vanishing values of ft or cr, the 
only possible zero in the cosine difference term (for the incompressible 
case) is at t =' 0. But at t = 0 a second-order zero appears in the 
numerator, completely overpowering the zero in the denominator. It 
should be noted that, as pointed out by Woolston and Runyan (ref. l) in 
the wall -interference problem, resonance is possible in the limit of 
zero Mach number if this is achieved by letting free-stream velocity go 
to zero while maintaining finite acoustic velocity. In the present case, 
however, incompressible flow in the classical sense is implied. 


It should be noted that, in the incompressible case. 



A 0 + 


J 0 (t) 


+ I|Aq + 


^Wt) + 


CXI 

He 

v=2 


v+1 



(B3) 


and F 0 (-k) is simply the same expression with -k inserted for t. 
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Now referring to KUssner and Schwarz (ref. ll) , or to Fung (ref. 12), 
it is found that in the single-wing case, for zero Mach number, under 
the trigonometric substitution of equations (4l) 

x = -cos 7 


^upper “ Slower 





00 

A n sin ny 

n=l 


If the upwash ratio 


v 

U 


is expressed in a cosine series 


y(y) 

U 


w. 


0 


w n cos n7 


n=l 


(B4) 


then there is a set of relations giving the pressure coefficients A n 
in terms of the (known) upwash coefficients as. follows: 



A 


1 


4 




Ag 

2 



- W2 


> (B5) 



ik 

n 




w. 


n 


(n ^ 2) 


where C = C(k) = F(k) + iG(k) is Theodorsen's well-known flutter func- 
tion. While equations (B5) constitute an effective solution in the usual 
sense, for the purposes of this investigation in combining effects with 
the cascade contribution, it is necessary to back track and find the 
contribution of each A n term to all the w n terms. Eventually coef- 
ficients of cos n 7 will be equated on both sides of equation (Bl) 

(which is equivalent to the Galerkin type orthogonalization of error 
with respect to cos my where m = 0, 1, 2, 3, • • •) and this requires 
a knowledge of the contribution of each component of pressure in 
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equations (4l) to every upwash component in equation (Bk) . The method 
of getting this contribution in the cascade portion is quite straight- 
forward and will be discussed at the conclusion of this appendix. The 
contribution of each pressure component to the downwash elements in the 
single-wing case requires a solution of the infinite system of equa- 
tions (B5) for the terms w n , assuming the terms An to be given. This 
development follows: 

First, it is noticed that the homogeneous system corresponding to 
equations (B5) is, for n ^ 2, 

0 = IT ( v n-l - v n+l) “ w n < b6 ) 

This is satisfied by Bessel functions Jn and Y n of the first and 
second kinds, respectively, with argument (-k) and coefficient (i)" n 
as may be seen by reference to the recurrence relations satisfied by 
the Bessel functions. Thus a general solution to equation (B6) is 

w n = ai -n J n (-k) + bi~ n Y n (-k) (B7) 

Now from equation (B7) a type of Green's function is constructed for 
the system of equations (B5) by imagining one of the terms A m to be 

unity while all others vanish. The corresponding terms w n are then 
precisely the desired elements of the single-wing upwash distribution 
corresponding to the term multiplied by A m in the pressure expansion. 
The first two equations of system (B 5 ) do not fall in the pattern of the 
Bessel function recurrence scheme and thus have the character of initial 
conditions in a difference equation. 

The construction of a Green's function proceeds by assuming first 

that 



in which case 
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a^ * 

J x [ik(C - 1) - c] - CkJ 2 


(B12b) 


the arguments (-k) being understood for the Bessel functions. In equa- 
tions (B12) the superscript ( 0 ) indicates that the values are for the 
particular case implied by equations (B 8 ) . Hence, equations (BIO) 
become, for this case, 

w n (°) = a^ 0 ) i“ n J n (-k) (n !> l) (B13) 


Note that Wq^ remains bounded as k -> 0, In fact, -* 


2C(o) 


1 

2 


as k-» 0 which checks the steady-state theory. Next, the term w n (-*-) 
is determined by letting 


Aq = a 2 = A 3 = \ = 


A 1 = 1 


0 


(BlA) 


in the system of equations (B 5 ) . All but the first two equations of 
equations (BIO) are satisfied if 

w n^ = a^ J n (-k)i" n (n ^ l) (B15) 


again dropping the Y n (-k) terms for reasons of boundedness of the 

terms ' w n ^ as k 0. The two remaining quantities a^ and Wq^ 

must be found so that the .first two equations of the system are satisfied. 
Thus 


0 = 


w, 


( 1 ) 


0 





ik 





(B1 6) 


and 
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“o' 11 - ^ - 1 ) 

aW . 111 

Jl(-k)^k - | + ij + ikJ 2 (-k) 


(BIT) 


which, together with equation (BI 5 ) completes this case. Once again, 

all w^ 1 ) terms are hounded as k-> 0. Continuing in this manner, 
that is, letting m run through the integers in the relation 

(m ± n) 

(m = n) 

it is found for m = 2 


A m ®mn ® 


A m = 5 !m " 1 


w n ( 2 > = i_n J n( -k ) 


w-^(2) 


w q (2) 


a( 2 ) kJ 2 (-k) 
k 


- k - i 


w x (2) 


1 

C 


a 


( 2 ) 


t (B18) 


(n ^ 2) 


y (B19) 


ik 2 J 2 (-k) 


2J 2 (-k) + kJ 5 (-k) + 


|-k-i 
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i“ n J n (-k) 


n 


(n £ 3) 


v- 


( 5 ) _ ik 2 a (3) J 3 (-k) 

k 2 + 2i/k - ]| + i 
0 


w (3) .affl A . i\ 

° ' 2 f °) 


V" = V J ' ( x - 5 + I) 


,(3) = v,(3) 
,(3) . 


51/2 


3J 5 (-k) + kJ^(-k) - 


(k - % + l)j 5 (-k) 


ik 2 ( k - + ij 

k 2 + 2i ^k - | + i 


For m = k, 

S h) =a(^) i- n J n (-k) 


(B20) 


w 


n 


(n £ h) 


v x (^) = 


k 5 a^ J^(-k) 


k 2 fk - | + i] + 31 


« 2 w . (1 - i + 1 


k 2 + 2ifk -|+i 


vJU - 
5 >2 


(4) 


k^ 

w^) 


w O (4) = 2 


k 2 + 2i/k - £ + i N 

l C > 


1 - i 


,( 4 ) _ 


v T.- ( -k ) - I 4 .J 1 f-k) + 

ik2 J 1+ (-k)[k2 + 2i\ 

(k - j 



— Ivu pr V *“■«. / M-O 2^ \ — -t*- / • 

k 2 | 

( k - 5 + 4 ) + 5i 

k 2 + 

2i^k 

-! +1 ) 

1 


(B21) 
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Finally, for all values of m 3> 2, a general expression is found for 
the solution to the system of equations (B5) with A n = 5nm* Letting 


a Q = 1 


M +1 ) 


a 2 = k2 aQ + 2ia^ 


a 3 = k^ a^ + 3ia-r 


“ta = k “ta-2 + 


gives the following general expressions: 


,<m) . 


(§)(D m - g 

mJ m (-k) + kj^jj-k) - ik 2 J m (-k) ^=2 

^-1 


(*) . a ( m ) (i)“ n J (-k) 


(n > m) 


(l ) _m j (_k) 


1 -i 


(m) ^1 


w 5 (m) 


(m) 


Wi, ( m ) 


= W-, (m) 22 

-L ~z 


w m -j (m) = w-i ( m ) 2m=Z 
ni— 1 -L -^in-2 
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Again, all are bounded as k-> 0. Thus, the contribution of 

each single-wing portion of the pressure series (eqs. (kl)) to every 
term in the upwash series (eq. (B4)) has been completely defined. 

There remains only to express, in correct form, the cascade contribu- 
tion. This is easily accomplished utilizing the relation (ref. 8 ) 

e iz cos 7 _ j o ( z ) _ 2jj 2 (z)cos 27 - J^(z)cos ky + . . . J + 

2i|j 1 (z)cos 7 - J^(z)cos 37 + . . . J 
00 

= ^ e v J v (z)cos vy (£24) 

v =0 

where e y = 1 when v = 0 and e v = 2(i) v when v > 0. Applying 

equation (B24) to equation (Bl) and using equations (B3), (B4), and 
(B 12 ) to (B23) gives 


.n =0 


with 


n=0 


V cos nr = V OOB nr 

Y W n (m) Am 


m=0 




Single-wing contribution 


jtk 

2 


00 

r — 1 


• ks 

sinh — 

) e n J n (k)cos n 7 

0 

1 

b - 1 

L—i 


cosh - cos (a + 

n =0 


b V b ) 


p N-k 

G 0 (t) - G 0 (-k) 

' -N-k 

t + k 


dT 


(B25) 


Gq(t) 


•i|T|F 0 (r) 


sinh 


rs 

b 


- 1 


cosh 

b 


cos a 


(0- 


TS ' 


Y e n J n ( -T ) cos 

n =0 


(B26a) 
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G n (-k) = - 


i|-k|F 0 (-k) 


cos a + 


- 1 \ e n J n (k)cos nr 


(B26b) 


The term F q (t) is given by equation (B3) and the term F 0 (-k) is 

given by the statement following equation (B3) . The Galerkin process 
now involves the retention of a finite number p of terms in equa- 
tions (4l) and the equating of coefficients of cos nr for the first 
p integers (n = 0, 1, 2, . . . (p - l) ) . This constitutes a system 
of p equations for p unknowns and may be carried as far as desired 
depending upon the available computing facilities. The major effort 
involved lies in the computation of the integral 


•N-k t + k 


(t) 

n' ' m' 


sinh 


cosh 


cos a - 




which, for a sizable number of parameter values, should be performed 
by a high-speed computer. 
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APPENDIX C 

COMPTJTATIOML PROCEDURES 


Following is a description of the computational procedure utilized 
in evaluating the cascade contribution to the integral equation ( 63 ) 
for purposes of solution by collocation. In order to compare the numer- 
ical results of the present method with those of the method of refer- 
ence 2, the special case of M = 0.5> zero stagger, antiphase motion 
(cr = 0, = l80°) , and gap-chord ratio of 3*80 was chosen. Calculations 
were made for eight values of k which are 0.02, 0.1, 0.2, 0.3> 0.4, 

0. 0.6, and 0.7- All computations were programmed for and executed 
on the Burroughs E-101 Electronic Digital Computer. 


The Burroughs Computer, as currently available at NYU, is of limited 
capacity having 100 12-digit memory locations and capable of accepting 
only 112 instructions which are introduced by means of pinboards . To 
multiply two numbers which are in different memory locations and write 
the result in a third location requires four instructions. The machine 
can accomplish approximately four multiplications a second and does not 
have floating decimal. Hence all problems must be carefully programmed 
and scaled prior to calculation. It is doubtful that this machine could 
be used for the general (nonzero -stagger and arbitrary -phase -motion) 
problem for, as will be shown, these two assumptions introduce several 
simplifications in the analysis. Even with these simplifications opti- 
mum use of the computer capabilities was required. 

In order to put the integral of equation ( 63 ) in a form more amena- 
ble for calculation the following adjustments were made: 


(l) To eliminate the necessity of programming square roots, sub- 
stitutions 


t = — cosh A 


t = — sin A ^ 


t = ■^5 cosh A 
P 2 


for 



zM < T < M, and T < d35 

P 2 = = P 2 P 2 ■ 


(Cl) 


were made, and 


respectively. 
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the range of integration -N ^ g t g N - — was subdivided into the 

p P 2 

three ranges indicated in equations (Cl). 

( 2 ) After introducing the factor exp ( x ) (see eq. (63)) into 

VP 2 / 

the integrand, the integral was separated into real and imaginary parts. 

(3) The factor G( — — | being unnecessary in the two ranges 

1 ^ -of-* 

- — < r -g — and — •g t -g /n — — the integral P — 

P 2 P 2 P 2 l P 2 / J 


was 


T + 4 


removed and evaluated analytically for these ranges. 

After carrying out these operations and introducing 


M^ = cosh' 


-1 


M2 = cosh 


-1 


eb -4 

kM \ p 2 J 

,2 / 




kM 


E 


-JO 
11 ?, 


= sin(-kx)F 0 /^ £ 

\p *} M 


sinh/l& 




- 1 


cosh 

b 


cos/^^ + 0\ 

\P 2b ) 


^l=r\ = cos(kx)F 0 ^\ | 
\P 2 / \P‘ / 


sinhj 

fsk' 
V b 


cosh - cos| 

^ + n\ 

^P 2 b ) 


- 1 


the integral in equation (63) was separated into the following parts: 
For range ^ < t < ( N - | , 

p 2 " “ \ p 2 / 


“ T Ml - Sln h2A - sin/™* cosh X + S^Vof- 
2 P J 0 cosh X + 1 \p 2 P 2 ) \P 2 


cosh A 


sinh/^^ sinh %) 

W L 


- 1 


cosh(|^ sinh A] - cos j cosh A - H 


/s*hM 

'Ip 2 * 


ax (02) 
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ikM 

2P 


X 


Mr 


sinh A 

cosh A - i 
M 


<>cos 


-kMx 


cosh A +. 


M-L±)f [d21 cosh A) sinh A 
P 2 ) > 2 j 


sinh I 

sinh a! 
(bp J 

\ 

sinh >' 
1 

— 

| - cosj 

(sfkM cosh x + fi \ 

j 


- 1 




> dA 


(C7) 


All integrals were evaluated using the Gaussian nine-point integra- 
tion (ref. 13, p. 285) and were checked by the Newton-Cotes closed- 
quadrature formula (ref. 13, p. 123) • The Gaussian method was accom- 
plished in all cases, except in the nearly resonant case, k = 0.7 and 

< A ^ 4., by subdividing the interval of integration into the nine 

prescribed points. For the Newton-Cotes method, the interval of inte- 
gration was divided into 30 equally spaced points and the six-interval 
formula (ref. 13, p. 123) was used five times. Because of the sharp 
peaking of the integrand in the above-mentioned case, the interval 

— < A < £ was divided into three segments and the Gaussian nine-point 

integration was carried out over each segment. This was checked again 
using the Newton-Cotes six-interval formula over 30 equally spaced points. 


In order to evaluate the integrand at each point it was necessary 
to generate sines, cosines, hyperbolic sines and cosines, and four Bessel 
functions, Jq, J^, Jg, and J^. The standard method of storing coef- 
ficients of the power series expansions, as recommended by Burroughs 
Corporation, would require about 70 of the 100 available memory locations 
in order that the truncated series used be accurate over the range of 
arguments of these functions. This was clearly impossible. A method 
was devised for generating the series by using a recursion formula for 
the coefficients that would serve for all the above series. For any of 
these series the nth term is obtainable from the (n - l) term by multi- 

plication by a factor of the form — -, where X is the argument! 

2n(2n + H) 

the plus sign is used for sinh X and cosh X and the minus sign, for 
all others; H = 1 for sin X and sinh X, -1 for cos X and cosh H, 
and 2v for J V (X) .• 
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The case considered of zero stagger (s* = 0 ) and antiphase motion 
(0 = ft = it) reduces the factor in the integrands of equations (C2) , (C 3 ), 
(C6), and (C 7 ) 


sinhf-^M sinh A 

Ibp 


coshh^^ sinh A 

W 1 


cos I?™, cosh A . a| 

W 2 / 


sinh/^l sinh a'’| 

L 1 (C8) 

coshf^^ sinh A | + 1 

\*P / 


Since the upper limit of integration was fixed by insuring that expres- 
sion (C8) be sufficiently small, accuracy would be lost by having to 
subtract the two very close numbers in equation (C8) when A is near 
the upper limit. However, this was overcome by using the identity 


sinh X 
cosh X + 1 


- 1 = 


cosh — / 

2 l 


-1 

y 

sinh — + cosh 




(C9) 


This identity not only permitted more accurate generation of sinh and 
cosh for smaller arguments ^ instead of X^ but it also retained 

accuracy in the evaluation of the term (C8) . 


For the general nonzero stagger and arbitrary phase condition, the 
problem will clearly not permit this elegant simplification and hence 
would exceed the capacity of the E -101 computer. 

In all calculations (series evaluations and infinite integral trun- 
cations) a maximum error of 0.0001 was permitted. 

In order to solve the integral equation (63) by collocation at 
three chordwise stations ^x = -i, 0, and i j, these three x-values were 
treated in the numerical evaluation of the integrals. 
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TABLE I 


TOTAL SINGLE-WING CONTRIBUTION TO REDUCED UFWASH V Q /U 
AT THREE CHORDWISE STATIONS FOR M = 0-5 


k 

X 

Aq Contribution 

Contribution 

A 2 Contribution 



Real 

Imaginary 

Real 

Imaginary 

Real 

Imaginary 

0.02 

0.5 

0.4519491 

0.0428256 

0.22578212 

0.02582204 

0.21671019 

-0.00488926 

.02 

0 

A515698 

.0487572 

.00906871 

.02734223 

.4332256 

- .0000756 

.02 

-.5 

.4511555 

.0546844 

-.20764554 

.02597331 

.21671191 

.00473800 

.10 

.5 

.5352984 

.1142825 

.26513120 

.0809435 

.21963380 

- .02613512 

.10 

0 

.5304671 

.1474325 

.04529458 

.0901383 

.4368039 

-.0018872 

.10 

-.5 

.5242158 

.1802806 

-.17464284 

.0847148 

.21983553 

.02236510 

.20 

.5 

.6450755 

.1297442 

.31652281 

.1172289 

.22520583 

- .056794409 

.20 

0 

.6338565 

.2054147 

.09028736 

.1397294 

.4447680 

-.0075312 

.20 

-.5 

.6160410 

.2794763 

-.13674993 

.1322414 

.22681011 

.041806889 

.30 

.5 

.7510416 

.0955288 

.36804027 

.1334792 

.23057684 

-.09213552 

.30 

0 

.7379910 

.2230042 

.1346790 

.1734030 

.4549473 

-.01687482 

.30 

-.5 

.7079702 

.3468238 

-.10136703 

.1669715 

.23595639 

.058772773 

Ao 

.5 

.8455129 

.0232180 

AI796165 

.1350070 

.23398967 

. .131978574 

AO 

0 

' .8394037 

.2106121 

.178177 

.1962812 

.4662342 

-.0298216 

AO 

-.5 

•7995035 

*3925083 

- .06791077 

.1938359 

.24663393 

.07355509 

• 50 

.5 

.9222588 

-.0798821 

.46486215 

.1245734 

.23402381 

-.17589716 

• 50 

0 

.9352460 

.1734976 

.2204926 

.2107855 

.4777981 

-.0462385 

• 50 

-.5 

.8901497 

A2128385 

- .03602700 

.2150847 

.25844802 

.08636572 

.60 

.5 

.9761701 

-.20748176 

.50753505 

.1041001 

.22949481 

-.22325336 

.60 

0 

1.0230507 

.1151938 

.26134806 

.2183600 

A889683 

-.0659546 

.60 

-.5 

.9794408 

.43592156 

-.00549879 

.2319856 

.27112929 

.09736951 

• 70 

.5 

1.0031925 

-35347302 

•54495248 

.0751580 

.21942893 

- .27323354 

• 70 

0 

1.1006352 

.0385219 

.30047885 

.2199862 

A891746 

-.0887732 

• 70 

-.5 

1.0669196 

A3825568 

.02381299 

.2453494 

.28448357 

.10669571 
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TABLE II 

EXTERNAL BLADE CONTRIBUTION TO UFWASH V Q /U AT 
THREE CHORDWISE STATIONS (TOTAL CONTRIBUTION 
LESS SINGLE-WING CONTRIBUTION) 

[m = 0.5; | = 7-60; s* = 0; a = l80^j 


k 

Y 

Aq Contribution 

Contribution 

A 2 Contribution 

(a) 

A 

Real 

Imaginary 

Real 

Imaginary 

Real 

Imaginary 

0.02 

0.5 

-0.05298 

-0.01785 

-0.01254 • 

-0.00904 

-0.004045 

0.000098 

.02 

0 

-.02514 

-.0185? 

-.00847 

-.00957 

-.004124 

•000070 

.02 

-.5 

-.01700 

-.01824 

-.00440 

-.00918 

-.004045 

.000041 

.1 

.5 

-.07794 

-.00621 

-.05512 

-.00449 

- .005992 

.001519 

.1 

0 

-.07052 

-.00899 

-.05108 

-.00574 

-.004141 

.001176 

.1 

-.5 

-.06219 

- .01148 

-.02695 

-.00685 

-.004122 

.001031 

.2 

.5 

- .09096 

.05060 

-.04511 

.01255 

- .002640 

.002695 

.2 

0 

- .O86O8 

.02504 

-.04018 

.00992 

-.005055 

.002515 

.2 

-.5 

-.08026 

.01985 

-.05697 

.00751 

-.005245 

.002507 

.5 

.5 

- .08515 

.05858 

-.04105 

.02680 

-.000911 

.002559 

.5 

0 

- .08215 

.05562 

-.05982 

.02424 

-.001556 

.002589 

.5 

-.5 

-.07974 

.04852 

-.05806 

.02167 

-.002021 

.002576 

a 

.5 

-.07748 

.07917 

-.05895 

.05854 

-.000242 

.001262 

.1+ 

0 

-.07796 

.07715 

-.05829 

.05708 

-.001056 

.001654 

a 

-.5 

. -.07669 

.07425 

-.05695 

.05529 

-.001695 

.001951 

.5 

.5 

-.08155 

.10885 

-.04006 

.05557 

-.001189 

-.000562 

.5 

0 

-.08025 

.iio86 

-.05842 

-05552 

-.002115 

.000270 

.5 

-.5 

-.07701 

.11119 

-.05597 

.05489 

-.002867 

.001021 

.6 

.5 

-.1025 

.18419 

-.04655 

.09567 

-.005269 

-.002861 

.6 

0 

-.09555 

.19156 

-.04079 

.09772 

-.006564 

-.001212 

.6 

-.5 

-.08175 

.19568 

-.05405 

.09815 

-.007205 

.000415 

• 7 

.5 

- .24585 

.74995 

-.08450 

.58755 

-.05865 

-.009206 

• 7 

0 

-.16857 

.77459 

-.04460 

.59526 

-.04055 

-.002671 

.7 

-.5 

- .08955 

.78649 

- .00405 

.59261 

- .04099 

.005944 


a Critical k for lowest resonance = 0. 715973* 
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TABLE III 

INFLUENCE COEFFICIENTS FOR LIFT AND MOMENT 
OF SINGLE WING AT M = 0-5 


[r indicates real; I indicates imaginary] 


k 


(0) A 0 

(0)*1. 

(0) A 2 

(1) A 0 . 

aA 

(l) A 2 

mm 

R 

2.1922 

0.00245 

0.00021 

-0.05923 

2.3109 

0.00000 

mm 

I 

-.2380 

.0604 

-.00002 

-.1188 

-.00055 

-.01538 

mgm 

R 

1.7846 

0.03141 

-0.00397 

-.2718 

2.3276 

0.00037 

■91 

I 

-.5248 

.2809 

-.OOO65 

**.2551 

-.01438 

-.07733 

0.2 

R 

1 .4902 

0.07702 

0.01188 

-.4285 

2.3565 

-0.00206 

.2 

I 

-.5737 

.5341 

-.00309 

-.2645 

-.04754 

-.1559 

0.3 

R 

1.3274 

0.1297 

0.02527 

-0.5218 

2.3901 

-0.00580 


I 

-.5 688 

.7801 

- .00744 

-.2421 

- .09112 

-.2358 

0.4 

R 

1.2258 

0.1953 

0.04357 

-0.5848 

2A280 

-0 .01212 

.4 

I 

-.5610 

1.0263 

-.01443 

-.2148 

-.1435 

-.3172 

0.5 

R 

1.1555 

0.2724 

0.06683 

-0.6327 

2.4703 

-0.02148 

.5 

1 

-.5604 

1.2751 

-.02487 

-.l8f4 

-.2054 

-.4005 

0.6 

R 

1.1015 

0.3716 

0.09514 

-0.6723 

2.5161 

-0.03458 

.6 

I 

-.5681 

1.5267 

-.04001 

-.1601 

-.2780 

-.4857 

0.7 

R 

1. 0511-2 

0.4967 

0.1347 

-0.7081 

2-5774 

-0.05492 

•7 

I 

-.5833 

1.7825 

-.0644 

-.1278 

1 

-.3648 

-.6017 



l/(*p 0 w 2 ) 


L 

[(0*0 + 

2 

H 

(1 - ika) 

O 

+ 

<? 

<3, 

+ ik 

U + 

1 


J 

L 

L c J 


L 2 jj 


M ea/(*Po' faSu2 ) = (a L /( It Po' bu2 )] 



(0) k 2 
2 J 


l 


+ e 


2 J 


(o)*o + 


f [(1)*© 


lira i 
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TABLE IV 

INFLUENCE COEFFICIENTS FOR LIFT? AND MOMENT OF BLADE 
ROW AT ZERO STAGGER, GAPCHOKD RATIO = 5.80, 

M » 0.50, AND o « fi = 180° 


[r indicates real; I indicates imaginary J 


k 


(o)^o 

(0) A 1 

(0)*S 

(l)Ao 

(D A 1 

(1)*2 

0.02 

R 

2.3328 

0.08788 

0.00217 

- 0.00957 

2.5555 

0.00012 

.02 

I 

-.1740 

.05666 

.00501 

-.08549 

-.00519 

-.01295 

0.1 

R 

2.0249 

0.1107 

«o. 00505 

-0.1692 

2.5674 

-0.00028 

.1 

I 

-.6359 

.2790 

-.0024 

-.5090 

-.01940 

-.07837 

0.2 

R 

1.6254 

0.1626 

0.01354 

- 0.3829 

2.4008 

- 0.00245 

.2 

I 

-.7718 

.514-28 

-.00639 

-.3630 

-.05152 

-.1577 

0.5 

R 

1.5952 

0.2524 

0.02799 

- 0.5164 

2.4472 

- 0.00629 

.5 

I 

-.7648 

.8086 

-.01529 

-.3410 

-.09347 

-.2592 

0.4 

R 

1.2637 

0.3^55 

0.04898 

-o. 6 o 84 

2.5167 

-0.01322 

.4 

I 

-.7649 

1.0831 

-.02571 

-.5159 

-.1526 

-.3240 

0.5 

R 

1.1670 

0.51459 

0.07645 

-0.6994 

2.6214 

-0.02621 

*5 

I 

-.8171 

1.3654 

- .04978 

-.3013 

-.2500 

-.4148 

0.6 

R 

1.0061 

1.0032 

0.1043 

-0.8538 

2.7894 

- 0.05746 

.6 

I 

-.9864 

1.5893 

-.1081 

-.2883 

-.4826 

-.5163 

0.7 

R 

-0.1822 

1.7376 

- 0.1246 

- 1.2242 

2.2158 

-0.2135 

.7 

I 

-.8514 

-.2052 

-.2255 

.2815 

-1.5488 

-5549 


l/(* Po UJ 2 ) = £ [ik (0) A 0 + 


+0(1- ika) 


[<o*> - 


+ ik 


(i) A o + ^ 




Mea/(«Po^ 2 U 2 ) = [aL/(* Po Wj2)] + | ^ (0)^) + • + 5 (0)*0 + M^J + • 
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Resonance lines 


Range in which ft is 
restricted to lie by- 
periodicity of problem 


(III 


Figure 2.- Relationship between resonance and singularities of cascade 
term. Letters are defined as follows: zone A, no poles of cascade 

term lie on real axis in t - plane ; zone B, t®, both lie on real 

axis in t - plane; zone C, rj*, both lie on real axis; zone D, 

T l’ T< 2* T ±’ T ? 1 ^ ° n rea ^ ax ^- s > an< 3- so forth; line aa, = t 

1 -1 

two poles coincide on real axis; line bb, T X = r 2 } ' fcwo P 0 "*" 63 c °i n “ 

cide on real axis; line cc, r" 1 = t*, two poles coincide on real 
axis; and so forth. 









‘feet of cascade or of • 
dtching about midchord 
lent amplitude in cascadi 

'lade ; a = l80°; M = 0.5 
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Figure 7 • - Effect of cascading or of tunnel walls on lift amplitude 
for pure translatory motion. JL^J = Lift amplitude in cascade; 

jL^'j = Lift amplitude of isolated blade; a = 180°; M = 0.5; 
s/b = 7 . 60 ; s* = 0 . 




lire 8.- Effect of cascading or of t 
for pure translatory motion. = 

= Lift phase angle of isolated 
s/b = 7.60: s* = 0. 







ure ±0.- Miect oi cascading or of tunnel walls on moment (about 
midchord) phase angle for pure translatory motion. = Cascad 

moment phase angle; 0^' = Single -wing moment phase angle; 0 = 1 

M = 0-5; s/b = 7.60; s* = 0. 





